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A recent paper [I. Klebanov et al. Mod. Phys. Lett. B 22 (2008) 3153; larXiv:0712.0433l 
claims that the exact solution of the Percus— Yevick (PY) integral equation for a system 
of hard spheres plus a step potential is obtained. The aim of this paper is to show that 
Klebanov et al.'s result is incompatible with the PY equation since it violates two known 
cases: the low-density limit and the hard-sphere limit. 
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1. Introduction 

Given a fluid of particles interacting via a certain potential 4>(r) , access to its (equi- 
librium) structural and thermodynamic properties is usually obtained by means of 
approximate integral equations,^ whose solution typically requires hard numerical 
work. Excep tion s are practically restricted to the Percus- Yevick (PY) equation for 
hard sphered^^ and sticky hard spheres,^ and the mean spherical approximation 
(MSA) for the hard-core Yukawa potential.^ 

The simplest potential including an energy scale and two length scales is that 
of hard spheres plus a step-function tail: 

{oo, r < a, 
e, a<r<\, (1) 
0, r>\. 

If e < 0, the step potential is attractive and Eq. (JXJ) describes the well-known square- 
well potential. On the other hand, e > defines the square-shoulder potential. The 
pure hard-sphere fluid is recovered if either e = (at arbitrary A/cr), or A = a (at 
arbitrary, but finite, e), or e — ► oo (again, at arbitrary A/cr). 

Starting from the PY integral equation, WertheinPl W as able to express the 
Laplace transform 

/•OO 

G(t)= / dre- tr rg(r), (2) 

J a 

where g(r) is the radial distribution function, in terms of quantities involving the 
cavity function y(r) = e^ r ^ kBT g{r) (where ks is the Boltzmann constant and T 
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is the temperature) only in the interval < r < A: 

_ (1 + ATtpK)t- 2 - F{t) + 2npt- 1 [Y{-t) - Y{t)\ 
[, ~ l + 2TT P t- 1 [F(-t)-F(t)} ' [> 

where p is the number density, K = —F'(0), and 

F(t) = - [ dre- tr rf(r)y(r), (4) 
Jo 

Y(t) = - dre~ tr dr' r' f(r')y(r')(r - r') [1 + f(r - r')]y(r - r'). (5) 

JO Ja+r 

Here, /(r) = e~^l kBT - 1 is the Mayer function. Equations Q-© apply not 
only to the potential {1} but more in general to any interaction with a hard-core at 
r = a and a finite range at r = A. 

In a recent paperf^ Klebanov et al. claim that they obtain the exact solution 
of the PY integral equation for the potential {T]). According to their approach, Eq. 
^ is complemented by 

y(r) = C x + C 2 r + C 4 r 3 , < r < A, (6) 

where the coefficients Ci, C2, and C4 are the solutions of a closed set of equations. 
Inserting Eq. ([6]) into Eqs. ([4]) and (O, one gets F(t) and Y(t), and hence G(t) 
through Eq. ©. 

The aim of this paper is to show that, in contrast to what is claimed in Ref. 
Eq. ([6]) is not compatible with the PY solution because it contradicts known results 
in the low-density limit as well as in the hard-sphere limit. 



2. Low-density limit 

The virial expansion of the cavity function is 

00 

y(r) = l + J2p n Vn(r), (7) 

71 = 1 

where the functions y„(r) are represented by sums of diagrams.^ In the special case 
of the potential ([I]), the first-order contribution yi(r) is given bj0 

yi (r) = (1 + 7 ) 2 ^, ff (r) - 2 7 (1 + 7)$<r,A(r) + 7 2 $a,a(0, (8) 

where 7 = e _e / fcsT — 1 and 

$ 0i6 (r) = -J- [3(a + 6) 2 - 2(b - a)r - r 2 ] (b - a - r) 2 9(6 - a - r) 

— — ^— [3(6 - a) 2 - 2(a + 6)r - r 2 l (a + 6 - r) 2 6(a + 6 - r), (9) 



I 
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Q(x) being Heaviside's step function. More explicitly, in the interval < r < 2a 
one has 

yi (r) = (1 + 7) 2 ^(4a + r)(2a - rf + 7 2 y^(4A + r)(2A - rf - 2 7 (1 + 7 ) — a 3 
s fo, 0<r<A-cr, , x 

+27(1+7) „, , ^2 „, , 2 _ _ . , 10 ) 

I 7r » ' — '- (A — a — r) , A — a < r < 2a. 

We see that, while yi(r) is a cubic function in the interval < r < A — a, it is 
a quartic polynomial function divided by r for r > A — cr. Moreover, the second 
derivative of yi(r) is discontinuous at r = A — a. Therefore, Eq. ^ is inconsistent 
with Eq. (|10p to first order in density. Since the PY theory yields the exact y\(r) 
for any interaction potential^ we conclude that Eq. ([6]) is inconsistent with the 
true PY solution. 



3. Hard-sphere limit 

As an independent test, let us now take the limit e — > keeping X/a fixed. In that 
case, as said above, the potential ([I]) reduces to that of hard spheres of diameter 
a, the value of A > a not play ing any role. The exact solution of the PY equation 
for hard spheres is well known JlEEE] j n particular, the functional form of y(r) for 
< r < 2a iP 

\A Q + A ir + A 3 r 3 , 0<r<a, 
y(r) = < (11) 
\^T e ' Kir + ^-e~ K2r cos (ujr + ip) , a<r<2a, 

where the coefficients A^, Bi, m, u), and (p are functions of density whose explicit 
expressions will not be needed here. It is quite clear that Eq. ([6]) cannot reduce to 
Eq. (fTTj) if e — > with A > a. Therefore, Eq. (J6j) is again incompatible with the PY 
equation. 



4. Conclusion 

In summary, the approximation presented in Ref. ^ is not the solution of the PY 
integral equation for the interaction potential (Q]), in contrast to what is claimed 
by Klebanov et al. The flaw in their derivation could be due to the fact that, at 
a given point, the authors discard the product G(t)F(—t) when manipulating Eq. 

Therefore, what they obtain is, at most, an approximation to the true solution 
of the PY equation, differing from the latter even in the low-density and in the 
hard-sphere limits. 
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